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fSJ , Abstract 

We consider the Dirac equation in the magnetic-solenoid field (the field of a 
. solenoid and a collinear uniform magnetic field). For the case of Aharonov-Bohm 

| solenoid, we construct self- adjoint extensions of the Dirac Hamiltonian using von 

7— I ' Neumann's theory of deficiency indices. We find self-adjoint extensions of the Dirac 

Hamiltonian and boundary conditions at the AB solenoid. Besides, for the first 
time, solutions of the Dirac equation in the magnetic-solenoid field with a finite 
CN \ radius solenoid were found. We study the structure of these solutions and their 

dependence on the behavior of the magnetic field inside the solenoid. Then we 
Q . exploit the latter solutions to specify boundary conditions for the magnetic-solenoid 

1 ' field with Aharonov-Bohm solenoid. 

<n : 
o . 

^ ; 1 Introduction 

PS. 

The present article is a natural continuation of the works [1, 2, 3, 4] where solutions of the 
Schrodinger, Klein-Gordon, and Dirac equations in the superposition of the Aharonov- 
Bohm (AB) field (the field of an infinitely long and innnitesimaliy thin solenoid) and a 
^ ■ collinear uniform magnetic field were studied. In what follows, we call the latter super- 

position the magnetic-solenoid field. In particular, in the paper [4] solutions of the Dirac 
equation in the magnetic-solenoid field in 2 + 1 and 3 + 1 dimensions were studied in 
detail. Then, in [5], these solutions were used to calculate various characteristics of the 
particle radiation in such a field. In fact, the AB effect in synchrotron radiation was 
investigated. However, a number of important and interesting aspects related to the rig- 
orous treatment of the solutions of the Dirac equation in the magnetic-solenoid field were 
not considered. One ought to say that in the work [4] it was pointed out that a critical 
subspace exists where the Hamiltonian of the problem is not self-adjoint. But the corre- 
sponding self-adjoint extensions of the Hamiltonian were not studied. The completeness 
of the solutions was not considered from this point of view as well. 

One has to remark that even for the pure AB field it was not simple to solve the 
two aforementioned problems. First, the construction of self-adjoint extensions of the 
nonrelativistic Hamiltonian in the AB field was studied in detail in [6]. In the work [6] 

*Dept. Fisica e Quimica, UNESP, Campus de Guaratingueta, Brazil; On leave from Tomsk State 
Pedagogical University, 634041 Russia; e-mail: gavrilovsp@hotmail.com 
^e-mail: gitman@dfn.if.usp.br 
■^e-mail: saa@dfn.if.usp.br 



1 



solutions in the regularized AB field were thoroughly considered as well. The need to 
consider self-adjoint extensions of the Dirac Hamiltonian in the pure AB field in 2 + 1 
dimensions was recognized in [7, 8]. The interaction between the magnetic momentum of 
a charged particle and the AB field essentially changes the behavior of the wave functions 
at the magnetic string [8, 9, 10]. It was shown that a one-parameter family of boundary 
conditions at the origin arises. Self-adjoint extensions of the Dirac Hamiltonian in 3 + 1 
dimensions were found in [11]. The works [12, 13] present an alternative method of treating 
the Hamiltonian extension problem in 2 + 1 and in 3 + 1 dimensions. It was shown in 
[14] that in 2 + 1 dimensions only two values of the extension parameter correspond 
to the presence of the point-like magnetic field at the origin. Thus, other values of the 
parameter correspond to additional contact interactions. One possible boundary condition 
was obtained in [9, 15, 16] by a specific regularization of the Dirac delta function, starting 
from a model in which the continuity of both components of the Dirac spinor is imposed 
at a finite radius, and then this radius is shrunk to zero. Other extensions in 2 + 1 
and 3 + 1 dimensions were constructed in the works [17, 18, 19] by imposing spectral 
boundary conditions of the Atiyah-Patodi-Singer type [20] (MIT boundary conditions) 
at a finite radius, and then the zero-radius limit is taken. In the works [21, 22] it was 
shown that, given certain relations between the extension parameters, it is possible to 
find the most general domain where the Hamiltonian and the helicity operator are self- 
adjoint. The bound state problem for particles with magnetic moment in the AB potential 
was considered in detail in the works [23, 24, 25]. The physically motivated boundary 
conditions for the particle scattering on the AB field and a Coulomb center was studied 
in [26]. 

The study of similar problems in the magnetic-solenoid field is a nontrivial task. In- 
deed, the presence of the uniform magnetic field changes the energy spectrum of the 
spinning particle from continuous to discrete. Thus, the boundary conditions that were 
obtained for a continuous spectrum cannot be automatically used for the discrete spec- 
trum. By analogy with the pure AB field it is important to consider the regularized 
magnetic-solenoid field (we call the regularized magnetic-solenoid field the superposition 
of a uniform magnetic field and the regularized AB field). Here one has to study solutions 
of the Dirac equation in such a field. The latter problem was not solved before, and is of 
particular interest regardless of the extension problem in the AB field. One ought to say 
that the Pauli equation in the magnetic-solenoid field was recently studied in [27, 28]. 

In the present article we consider the Dirac equation in the general magnetic-solenoid 
field (the uniform magnetic field and the AB field may have both the same and opposite 
directions) and in the regularized magnetic-solenoid field. First we construct self-adjoint 
extensions of the Dirac Hamiltonian using von Neumann's theory of deficiency indices. We 
demonstrate how to reduce the (3 + 1) -dimensional problem to the (2 + l)-dimensional 
one by a proper choice of the spin operator. We find self-adjoint extensions of the Dirac 
Hamiltonian in both above dimensions and boundary conditions at the AB solenoid. 
Then, we study properties of the corresponding solutions and energy spectra. We discuss 
the spectrum dependence upon the extension parameter. In the regularized magnetic- 
solenoid field, we find for the first time solutions of the Dirac equation. We study the 
structure of these solutions and their dependence on the behavior of the magnetic field 
inside the solenoid. Then we use these solutions to specify boundary conditions for the 
singular magnetic-solenoid field. To this end, we consider the zero-radius limit of the 
solenoid. One ought to say that the problem of the Hamiltonian extension in a particular 
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case of the magnetic-solenoid field (both fields have the same directions) was considered 
in [29] (scalar case) and in [30] (spinning case in 2 + 1 dimensions). However, the 3 + 1 
dimensional spinning problem was not studied as well as the relation of the extensions 
with the regularized problem. 



2 Exact solutions 

Consider the Dirac equation (c = h = 1) in (3 + 1) and (2 + 1) dimensions, 

id ty = HV, H = 7° ( 7 P + M) . (1) 

Here Y = (7°,7) , 7 = (7 fc ), P k = id k - qA k , k = 1,2, for 2 + 1, and k — 1,2,3, for 

3 + 1, v — (0, k) ; q is an algebraic charge, for electrons q = — e < 0. As an external 
electromagnetic field we take the magnetic-solenoid field. The magnetic-solenoid field is 
a collinear superposition of a constant uniform magnetic field B and the Aharonov-Bohm 
field B AB (the AB field is a field of an infinitely long and infinitesimally thin solenoid). 
The complete Maxwell tensor has the form: 

F Xu = B(6l6l-6l6l),B = B AB + B. 

The AB field is singular at r = 0, 

B AB = ^(x^ix 2 ) . 

The AB field creates the magnetic flux <3>. It is convenient to present this flux as: 

$ = (Z + fi) $0, $0 = 27r/e , (2) 

where l is integer, and < /i < 1. 

If we use the cylindric coordinates if, r : x l = r cos (p, x 2 = r sin ip, then the potentials 
have the form, 

a a n a 1 M S ^ nL P a r, . , >, COSW 

A = 0, eA 1 = [l + fi + A r — ^, eA 2 — — [l + /i + A r ^ , 

r r 

(A 3 = in 3 + 1), A (r) = eBr 2 /2 . (3) 



2.1 Solutions in 2+1 dimensions 

First, we consider the problem in 2 + 1 dimensions. In 2 + 1 dimensions there are two 
non-equivalent representations for 7-matrices: 

7 ° = a 3 , 7 1 = ta 2 , 7 2 = -^C, C = ±1 , 

where the "polarizations" ( = ±1 correspond to "spin up" and "spin down", respectively, 
cr = (a 1 ) are Pauli matrices. In the stationary case, we may select the following form for 
the spinors ^f(x), 

= exp {-iex } tpP (x ± ) , C = ±1, x ± = (O,^ 1 ,^ 2 ) . (4) 
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Then the Dirac equation in both representations implies: 

(<rP ± + Ma 3 ) ^\x ± ) = e^\x ± ), P± = (0,P 1 ,P 2 ) , (5) 

(aVPid 1 + Ma 3 ) ^~ 1) (x±) = e^- 1] (x ± ) . (6) 

We remark that the energy eigenvalues can be positive, e = + e > , or negative, e = 
-E < 0. One can see that 

4- 1] M^a 2 ^l(x ± ). (7) 

Further, we are going to use the representation defined by ( = 1. 

As the total angular momentum operator we select J = —id v + a 3 /2 that is the 
dimensional reduction of the operator J 3 in 3 + 1 dimensions. The operator J commutes 
with the Hamiltonian H. Then the spinors ipP have to satisfy Eq. (5) and the equation 



j-j --)^w ( x±) , i eZ . 



(8) 



Presenting the spinors in the form 



4 1] (x±) = gi{ip)i>i (r) , gi(<p) = 7y= ex P (l - k ~ \ (1 + } , 



we find that the radial spinor ipi (r) obeys the equation 

hipi(r) = eipi(r), h = U + a 3 M, 



n = -i { d r + 



p + l--{l-a 3 )+A(r) 



(9) 

(10) 
(11) 



Here h is the radial Hamiltonian, n defines the action of the spin projection operator on 
the radial spinor in the subspace with a given /, 

<?P ±9i(<p)fpi (r) = gi(ip)Thl)i (r) . 

It is convenient to present the radial spinor in the following form 

^(r)= [a 3 (e-U) + M] Ul (r), (12) 

where 



M r ) = ^2 c ° u iA r ) > u iA r ) = <t>iA r ) v <>i 

<T=±1 
1 



Vi 





o r = I i 



(13) 



and c a are some constants. It follows from (10) that U 2 u = (e 2 — M 2 ) u, therefore the 
radial functions 4>i,a{ r ) satisfy the following equation: 



(f_ d__P 1 
P df? + Tp~A + 2 



Tp 



M0 = 0, (14) 



p = 7 r 2 /2, 7 = 6(51, £ = sgmB, u = p + I - (1 + a) /2, u = e 2 - M 2 
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Solutions of the equation (14) were studied in [4]. Taking into account these results, we 
get: 

For any /, there exist a set of regular 1 at r = solutions Jj(T = {4> m ,i,a- , m — 0, 1, 2, . . .), 

<Pm,lA r ) = Im+\v\,m(p) • (15) 

Here I n>m (p) are the Laguerre functions that are presented in Appendix A. 

For / = there exist solutions irregular at r = 0. A general irregular solution for / = 0, 
/i^O reads: 

a = »-(l + a)/2, 2A = w/ 7 -£[a*-(1-<t)/2] , (16) 

where W\ :a / 2 are the Whittaker functions (see [31] , 9.220.4). The spinors in (10) con- 
structed with the help of the latter functions are square integrable for arbitrary complex 
A. The functions ip\ tOC were studied in detail in [4], some important relations for these 
functions are presented in Appendix A. We see that interpretation of uj as energy is impos- 
sible for complex A. For real A there exist a set of solutions (16) which can be expressed 
in terms of the Laguerre functions with integer indices: 

C,+l( r ) = Im+ii-l,m(p) , <T = +1, 771 = 0, 1,2, 

C,-lW = Im-»,m{p) , <r=-l, m = 0,1,2,... . (17) 

All the corresponding solutions ipi(r) of Eq. (10) are square integrable on the half-line 
with the measure rdr. The Laguerre functions in Eqs. (15), (17) are expressed via the 
Laguerre polynomials. 

Eigenvalues uj and the form of spinors depend on sgn£>. Below we present the results 
for B > 0. The results for B < can not be obtained trivially from ones for B > 0. We 
present them in Appendix B. The spectrum of uj corresponding to the functions 4> m ,i,aij') 

TGctds 

J 2 7 (m + Z + /x), Z-(l + <r)/2>0 , , 

\ 2 7 (m + (l + <7)/2), Z-(l + <r)/2<0 ' K > 

and the spectrum of u corresponding to the functions 4>m,a( r ) reads 

u= | 2 7 (m + „) a = l 
[ 2 7 m, (j = — 1 

We demand the spinors ui (r) to be eigenvector for LT, such that the functions u m ^^± 
have to obey the equation 

n« m ,(,±W = ±v^m,±W • (20) 

Then we can specify the coefficients in (13). 
In the case uj — 0, 

""•' W = ( 1,,,-M ) ' ' S " l! ^»-(<..(D )•'-»• (21) 



1 Here we use the terms "regular", "irregular" at r = in the following sense. We call a function to 
be regular if it behaves as r c at r — with c > 0, and irregular if c < 0. We call a spinor to be regular 
when all its components are regular, and irregular when at least one of its components is irregular. 
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That can be easily seen from the relations (88) - (91) for the Laguerre functions I n ,m{p) . 
In the case uj ^ 0, 

W) = {^Mr))' l - hul = 2l(m + l+ " ) ' 
«™+i,i,±(r) = (^;W iW ), 1 <-i,.=2 7(m+ i), 

uUJr) = (^ i(r| ),H..W» + l), 

= (^.(r))-'- -"-^")- (22) 

For u; 0, we construct solutions of the Dirac equation using the spinors u correspond- 
ing to the positive eigenvalues of the operator II. These solutions have the form, 

ipm,i( r ) = N K {e-V^) + M] u mil ,+ (r), 1^0, 

lti'(r) = N [a 3 (e-V^)+ M] u££(r), 1 = 0, (23) 

where N is a normalization constant. Substituting (23) into (10) we obtain two types of 
states corresponding to particles + ip and antiparticles _i/j with e = ±e = ±V M 2 + uj, 
respectively. The particle and antiparticle spectra are symmetric, that is = for 
given quantum numbers m, I. 

Consider the case uj = 0. As it follows from (10), (21) only negative energy solutions 
(antiparticles) are possible. They coincide with the corresponding spinors u up to a 
normalization constant 

_^o,j(r) = Nu ,i(r), I < -1; _^(r) = Nu^r), 1 = 0. (24) 

Thus, only antiparticles have the rest energy level. The particle lowest energy level for 
I < is + e = ^M 2 + 27. 

All the radial spinors ipm,i{ r ) are orthogonal for different m. The same is true both 
for the spinors and ipm- m the general case, the spinors of the different types are not 
orthogonal. By the help of Eq. (94) of Appendix A, one can prove this fact and at the 
same time calculate the normalization factor which has the same form for all types of the 
spinors, 

N = - l = jfi . (25) 

\e - v^J) 2 + M 2 

Besides, on the subspace I = there are solutions of Eq. (10) that are expressed via 
the functions ip\, a (p) (16). We present these solutions as follows, 

ij> u (r) = [a 3 (e - n) + M] u„(r) , 

u u {r) = ciu W)+ i(r) + c_itt w _i(r), u^r) = 4>u,a{ r ) v ° ■ (26) 
Using the relations (100) for the functions -?/>A,a(p), we obtain the useful expressions 

n« W| i(r) = iy/Tr/Uu-^r), n Vi( r ) = -^-S=^,+i( r ) , ( 27 ) 
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By the help of Eq. (102) from Appendix A, one can see that the spinors ip u (r), Vv( r )> 
uj 7^ uj' are not orthogonal in the general case. 

The fact that on the subspace / = (in what follows we call this subspace the critical 
subspace, and the subspace / ^ the noncritical subspace) there exist solutions with 
complex eigenvalues indicates that the radial Hamiltonian is not self-adjoint, at least on 
this subspace. 



2.2 Solutions in 3+1 dimensions 



To exploit the symmetry of the problem under z translations, we use the following repre- 
sentations for 7- matrices (see [15]), 



7 



<j 
-a 3 



, 7 



ia 2 
-ta 2 



, 7 



-la 






ia 1 



, 7^ 



/ 
-/ 



In 3 + 1 dimensions a complete set of commuting operators can be chosen as follows 

(7 5 = -itVtV), 



H, P 3 = -id 3 , J 3 = -id v + S 3 /2, S 3 = 7 5 7 3 [M + 7 3 P 3 ) /M . 
Then we demand the wave function to be eigenvector for these operators, 

S 3 ^ = sM/Mm. 



(28) 



(29) 
(30) 
(31) 
(32) 



Here M = a/M 2 + (p^) 2 , p 3 is z-component of the momentum and j 3 is ^-component 
of the total angular momentum. Remark that the energy eigenvalues can be positive, 
e = + e > , or negative, e = s < 0. The eigenvalues j 3 are half-integer, it is convenient 
to use the representation: j 3 = (l — Iq — |), where / = 0, ±1, ±2, ... . To specify the spin 
degree of freedom we select the operator S 3 which is the z-component of the polarization 
pseudovector [35], 



(33) 



eigenvalues of the corresponding spin projections are sM/M, s = ±1. 

Then in 3 + 1 dimensions one can separate the spin and coordinate variables and get 
the following representation for the spinors 



exp {— iex° + ip 3 x 3 } ^ s (x±) 



N 



1 + [p 3 + sMj /Mj i> e , s (x ± ) 
-1 + (p 3 + sM) /M] i/> e , s (x ± ) 



(34) 



Here ip £ ,s( x ±) are two-component spinors, x± = (0, x 1 , x 2 , 0), N is a normalization factor. 
As a result, the equation (29) is reduced to the equation 



(aP ± + sMa 3 ) 4> £>s (x ± ) = e^ £>s (x ± ), P ± = (0, P u P 2 , 0) . 



(35) 
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Presenting %p £yS (x±) in the form 



(36) 



where g% (ip) is given by Eq. (9), one comes to the radial equation 



h a ipi, a (r) = eipi >9 (r), h s = U + sMa 3 , 



(37) 



where h s is the radial Hamiltonian acting on the subspace with the spin quantum number 
s, II is given by Eq. (11). We remark that 



One can see that at fixed s and p 3 , Eq. (35) is similar to Eq. (5) in 2 + 1 dimensions. 
Thus, after separation the angular variable with the help of (9), the radial spinor (36), 
V'J.+i ( r )) can be obtained from the radial spinor (9), tpi (r), with the substitution M by 
M. The same is true for the particular case / = 0. Here the radial spinor ip W) +i (r), can 
be obtained from the radial spinor (26), ip^ (r). 

Using the results for (2 + l)-dimensional case, one concludes that in the critical sub- 
space complex eigenvalues of Eq. (29) exist. That means the Hamiltonian in 3 + 1 
dimensions is not self-adjoint. 



As well-known [7, 8, 15], the radial Hamiltonian in the pure AB field requires a self-adjoint 
extension for the critical subspace 1 = 0. As a result [8] one gets a one-parameter family 
of acceptable boundary conditions. In the case of our interest, the external background 
is more complicated, it includes besides the AB field a uniform magnetic field. The wave 
functions and the spectrum in such a background differ in a nontrivial manner from ones 
in the pure AB field. Thus, the problem of self-adjoint extension of the Dirac Hamiltonian 
in such a background, which is considered below, is not trivial. 

3.1 Extensions in 2+1 dimensions 

First, we study the (2 + l)-dimensional case. To this end we use the standard theory of von 
Neumann deficiency indices [36]. The (2 + l)-dimensional case was formally considered in 
[30]. We reproduce calculations for the (2 + l)-dimensional case in terms of the functions 
of Sect. 2. We generalize the results of [30] for arbitrary sign of B that allows to determine 
the non-trivial spectrum dependence on the signs of B, $. The (2 + l)-dimensional case 
results are necessary in order to extend this result to the (3 + l)-dimensional case. 

The Hamiltonian (1) acts on the space of two-spinors ip (x±). The angular variable can 
be separated, Eq. (10), that enables us to single out the radial Hamiltonian h (10) which 
acts on the space of two-spinors ip (r) . We start with the choice of the domain of definition 
of h, T> (h). Let T> (h) be the space of absolutely continuous square integrable on the half- 
line (with the measure rdr) and regular at the origin functions. One can make sure that h 
is symmetric on the domain T> (h). To determine whether the Hamiltonian is self-adjoint 
we have to define its deficiency indices, n± (h) = dim (V^, T> ± =Kei{h) =)= iM), where W 



ipe,-i(x±) = crV- £ ,i(:r_L) . 



(38) 



3 Self-adjoint extensions 
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is the adjoint of h. That is we have to find the number of linearly independent solutions 
of the equations, 



/ity^r) = ±iMiP ± (r), h ] = n f + a 3 M , (39) 
tf = -iS^d r + y fx + l-^(l-a 3 )+A(r) | a 1 . (40) 

Here M is introduced by dimensional reasons. For both cases (r) ) there exist only 
one linear independent square-integrable solution, for I = 0, that reads, 



^(r) = N[ V; ± ; /4 ^, , , ),B>0, (41) 



±e ±-/4y^ _ i(r) 
0i(r) 



(0 = ± e ±W4^ _ l(r) J , B < , (42) 

where, using (16), 

Mr) = il>xM, 2A = "2M 2 /7 - £ Ox - (1 - a) /2) , <r = ±1 . 

Thus, on the non-critical subspace the deficiency indices are (0,0), and on the critical 
subspace the deficiency indices are (1,1). Therefore, on the non-critical subspace the 
radial Hamiltonian is (essentially) self-adjoint, and on the critical subspace the radial 
Hamiltonian has self-adjoint extensions. Besides, there exist the isometry from T> + into 
T>~ , il) + (r) — > e in ip~(r), (leK. According to von Neumann's theory, the extensions of 
a closed symmetric operator 2 are in one-to-one correspondence with a set of isometries. 
Thus, self-adjoint extensions of the original operator h form the one-parameter family 
labelled by the parameter fi, h n . The domain of h n reads, 

V (h n ) = { X (r) = i) (r) + c (^+ (r) + e in ^~ (r)) : ip (r) G V (h)} , c G C , (43) 

where x is a two component spinor, x = (x\x 2 )- The behavior of the functions from 
V (h n ) at r — > is defined by the behavior of x(r). Using the behavior (101) of the 
function ip\ t a(p) at small p, we find, 

xM^W 1 i = I ( tan o_ l)r0i)r(1+Ma/7) {, ) > ^ > u , 

{ (tan§-l)l»r(l^+M 2 / 7 ) V T / ' 

One can verify that the limit 7 — > of the right hand sides of (44) coincides with the 
corresponding expression obtained in [8] in the case of pure AB field. For our purposes it 
is convenient to pass from the parametrization by Q to the parametrization by the angle 
6, < < 2tt, such that 

lim / „A U = % tan - + — . 45 

r^o x 2 (r) (Mr)" V 4 2 / 



2 We remark, that every symmetric operator has a closure, and the operator and its closure have the 
same closed extensions [36]. 
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To guarantee the self-adjointness of the Hamiltonian one has to demand the functions of 
its domain to satisfy Eq. (45). 

Thus, the solutions (26) obtained in Sect. 2 must be subjected to the condition (45) 
at r -> 0. Then, with the help of (27), (109) and (101), we find 

tan ( * + ® ) = " M WMr(1 -^» V ; J ' ^ > U . (46 ) 
I 4 2 1 ) m r(i- M )r(i-W2 7 ) R n v ' 

\ (e-M) 2M-ir( M )r(l- M -^/2 7 ) ^ 7 y ' ° ^ U 

3.2 Extensions in 3+1 dimensions 

Now we pass to the (3 + l)-dimensional case. Usually, the helicity operator Sh = SP/|P| 
is used as the spin operator. It relates to t-component of polarization pseudovector (33) 
as, Sh = S°M/ |P|. In this case it is necessary to find a common domain for two operators: 
H and Sh- That is not a trivial problem even in the special case p 3 — [21, 22]. Moreover, 
not for all the extension parameter values of the Hamiltonian there exists a self-adjoint 
extension of the operator Sh- This is the principal reason of our choice S 3 as the spin 
operator. 

In (3 + l)-dimensional case the Hamiltonian (1) acts on the space of 4-spinors of the 
form (34). The Hilbert space of 4-spinors (34) can be presented as the direct sum of two 
orthogonal subspaces with respect to the value of the spin quantum number s: T> (H) = 
© Then we consider the Hamiltonian (1) on each of the subspaces. Using 

Eqs. (35), (9) allows to single out the radial Hamiltonian h s acting on the subspace with 
the spin quantum number s. 

We choose the domain of definition of h s , T> (h s ) as follows. Let T> (h s ) be the space 
of absolutely continuous square integrable on the half-line (with the measure rdr) and 
regular at the origin functions. The radial Hamiltonian h s is symmetric on the domain 
V (h s ). Now we apply von Neumann's theory of deficiency indices to each of the subspaces. 

To define the deficiency indices of operators h s we have to solve the problem, 

hl^jf (r) = ±isMipf (r) , h\ = n f + sMa 3 , s = ±1 , (47) 

where h\ is the adjoint of h s , W is given by (40). One can see that Eq. (47) is similar to 
Eq. (39). Then, using Eqs. (41) [or (42)], (38) one obtains that for I — the solutions 
read, 

^M = Jv(^ 4 ^_ i(r) ),B>0, (48) 
^M = ^(t + e-%0,,_ l( r))- B<O - < 49 > 

<Ps,o{ r ) = i>\a(p) , OL = V ~ (1 + O") /2 , 

2A = -2M 2 /7 - f (a* - (1 - <r) /2) , <t = ±1 . 

whereas for / ^ there exist no square integrable solutions. Therefore, for each subspace 
s = ±1 on the non-critical subspace the deficiency indices are (0,0), and on the critical 
subspace the deficiency indices are (1,1). Thus, on the non-critical subspace the radial 
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Hamiltonian h s is (essentially) self-adjoint, and on the critical subspace it has self-adjoint 
extensions. 

Using the results of (2 + l)-dimensional case we conclude that on each subspace s = ±1 
self-adjoint extensions of the radial Hamiltonian h s form the one-parameter family labelled 
by the parameter Vl S) h^ s . The domain of h^ 3 on each subspace reads: 

V (^ s ) = {Xs (r) = (r) + c (V+ (r) + e^J (r)) : ^ (r) G V %)} , c G C , (50) 

Using the parametrization by the angle Q s similar to (45) we define the condition for the 
functions from domain V (h^ s ) at r — > as follows 

lim v J ^ =zstan( T + — M, s = ±l. (51) 



(r) (Mr)" 



r-o 2 ^ /^JT V 4 2 



Therefore, in each subspace s = ±1 solutions ^ WjS (r) on the critical subspace must be 
subjected to the condition (51) at r — > 0. Thus, in 3 + 1 dimensions there exist the 
two-parameter family of self-adjoint extensions of the Hamiltonian. 

3.3 Spectra of self-adjoint extensions 

Let us study spectra of the self-adjoint extensions h n . To this end we have to solve the 
transcendental equations (46) for uj considering two branches of e, one for particles and 
another one for antiparticles, ±e = ±V M 2 + uj. Introducing the notations 

uj = 2 7 x, x = ,x = ( ? e 2 -M 2 )/2 7 , Q{x) = ^ + 1, ? = ±, 
we can rewrite Eq. (46) for £? > as follows, 

Having a; for B > 0, one can obtain for i? < making the transformation 

Therefore, below we consider the case B > only. 

Possible solutions x = x (rj) of the equation (53) are functions of the parameter 77 (of 
/i, 7/M 2 , 6) and are labelled by m = 0, 1, ... . One can find the following asymptotic 
representations for these solutions at — > , 

sin (7r/x) r(m + 1 — /i) 

x m (77) = m + Ax m , Ax m = — 77, m = 1,2,3,..., 

7iT(m)(y (to) 
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All x m (0) , m — 1, 2, ... are positive and integer. The asymptotic representation of + x (rj) 
at — > is discussed below. The function + x (rj) vanishes at the point 77 = 2T(fi) and, 
in the neighborhood of the latter point, has the form 

I \ r(/i) -?^/2 

+I »<'» = ro.)( W (55) 

Here t/>(x) is the logarithmic derivative of the gamma function T(x), and — 0.577 

is the Euler-Mascheroni constant [37]. At |??| — > oo we found the following asymptotic 
representations, 

<Pm (qr)) =m + fi + Ax m , m = 0, 1, 2, . . . , 77 -> oo , 

(a?) = m - 1 + A 4 + Aa; m , m = 1, 2, 3, . . . , i] -> -oo , 
sin (7r/x) r(m + /i)Q (m + /i) 
nF(m + 1)?7 

These approximations hold true only for |Ax m | <C A* , |x (r/)| <C /i. 

According to [38] (see there Theorem 8.19, Corollary 1) if 7\ and T 2 are two self- 
adjoint extensions of the same symmetric operator with equal finite defect indices (d, d) 
then any interval (a, b) C 1Z not intersecting the spectrum of Ti contains only iso- 
lated eigenvalues of the operator T 2 with total multiplicity at most d. Let us select 
the extension h n at 6 = n/2 with the eigenvalues + e = My/l + 27 + x (00) /M 2 and 
± £ = ±M a/1 + 27-ta; m (±00) /M 2 , m > 1. Then the above theorem implies that if (a, b) 
is an open interval where a, b are two subsequent eigenvalues of h n at = tt/2, or ±£ = 0, 
then any self-adjoint extension h at 7^ 7r/2 has at most one eigenvalue in (a,b). Ac- 
cording to [39] (see there Chapter VIII Sect. 105 Theorem 3) for any e G (a, b) there 
exist a self-adjoint extension /i n with the eigenvalue e. As it follows from (53), (56), on 
the ranges (m — 1 + fi < ±x m (77) < m + /i, m > 1), (— M 2 /2^y < + x (77) < /i) the 
functions ±x (rj) = (±s 2 — M 2 ) /2 / ~f are one-valued and continuous. This observation is in 
complete agreement with the above general Theorems. The functions ±x m (rj) were found 
numerically in the weak field, 7/M 2 <C 1, for some first m's. The plots of these functions 
(for \x = 0.8) see on Figs. 1 and 2. 

One can see that 5x m = x m+ i (rj) — x m (77) — > 1 with increasing m. It follows from the 
equation (53) that 

Sx m - 1 = 7T _1 {cot (nx m ) - cot [ir (x m - /i)]} -1 ( - — - - 5Q J , m > 1 , (57) 

where SQ = In | < l/x m . The curve x 5 (77) may give an idea how the functions 
a^m (?7) behave at big m . 

Below we discuss some limiting cases. 

Consider weak fields B, for which 7/M 2 <C 1, and nonrelativistic electron energies, 
%m (v) 7/M 2 <C 1. Here the functions ±x (77) change significantly in the neighborhood of 
77 = only. Beyond the neighborhood of 77 = the functions ±x (77) take the values close 
to the corresponding asymptotic values given by (56). 

In the ultrarelativistic (77) 7/M 2 1, the behavior of x m (rj) qualitatively 

depends on /i. One can distinguish three cases: /j, < 1/2, /j, > 1/2, /j, = 1/2. If /j, < 1/2 
then the interval near 77 = on which the functions change significantly diminishes with 
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Figure 1: Particle lowest energy levels in dependence on the parameter r] + = 
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Figure 2: Antiparticle lowest energy levels in dependence on the parameter r]^ = 

r(l-u) I 2M 2 J lcU1 U ^ 2 J 



m increasing. If ji > 1/2 then this interval grows with m increasing. For /i — 1/2 and 
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~~ 1 < (l + H) < I' we ^ ne as y m Ptotic representation, 

fa) = m + q + , m > 1 . (58) 

One can see that negative ±x fa) exist only for r\ > . That is, in the problem 
under consideration for rt/2 < O < 3tc/2 there exist only one particle state and only one 
antiparticle state with energies |e| < M. The same situation was observed in the pure 
AB field case [8]. The minimal admissible negative x fa) is defined by the condition 
e = 0. In strong fields B, for which 7/M 2 ~ 1, the quantity Xq fa) is close to zero. Let 
©o correspond to such an extension that admits e — 0. The value of 60 is defined by the 
expression 

t e,A r(i - rtrfc + A/727) /m 2 \'-" 

In weak fields, 7/M 2 C 1, i fa) take big absolute values, and the angle 6 is defined by 
the expression 

(tx 9 \ r(l-u) / 

and does not depend on the magnetic field. It follows from (59) that in the superstrong 
fields B, for which 7/M 2 3> 1, the angle 6 does not depend on the magnetic field as 
well. 

In weak magnetic fields, 7/M 2 <^ 1, and for nonrelativistic energy values, xq^/M 2 <^ 1, 
one can get relations 

+ x ( V ) = -(2/ V f^ , (61) 
_x ( V ) = -{ V M 2 / 1 ) 1/ > 1 (62) 

that are valid when 77 is small in (61) and r/M 2 /7 ^> 1 in (62). 

Let us consider the particular case = —n/2. It follows from (46) that for B > 0, 
there exists = — M. The energies \e\ > M are defined by poles of T(l — x) or of 
r(l — n — x) for B > or B < 0, respectively. The spectrum e coincides with one defined 
by Eqs. (107), (23) for ip 1 . Moreover, using the relation (98), we can see that the spinors 
ipujij-) coincide with ip 1 up to a normalization constant, 

i> u (r) oc ip 1 (r) for = -tt/2 . (63) 

In the case = n/2 we have the following picture: It follows from (46) that for B < 
there exists + e = M. The energies \e\ > M are defined by poles of r(/j — x) or of T(l — x) 
for B > or B < 0, respectively. The spectrum e coincides with one found by Eqs. 
(107), (23) for ip 11 . From (98) it follows that the spinor ip^r) coincides with ip n up to a 
normalization constant, 

ip u (r) oc i\) n (r) for = tt/2 . (64) 

Using results for B < which are presented in Appendix B one can conclude that the 
spectrum asymmetry takes place for the spinning particles in the magnetic-solenoid field. 
There is a relation between the three-dimensional chiral anomaly and fermion zero modes 
in a uniform magnetic field [32] (for review see [33, 34]). One can see the effect also takes 
place in the AB potential presence. 
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The spectrum asymmetry is known in 2+1 QED for the uniform magnetic field. In the 
uniform magnetic field the states with uj = for / 7^ are observed if sgn/ = — sgn£? (for 
antiparticle if B > and for particle if B < 0). The spectrum changes mirror-like with the 
change of the magnetic field sign. One can see that for Z 7^ the spectrum properties in 
the magnetic-solenoid field is similar to the spectrum properties in the uniform magnetic 
field. The presence of the AB potential is especially essential for the states with I = 0, 
when the particle penetrates the solenoid. 

Spectra in 3 + 1 dimensions can be obtained from the results in 2 + 1 dimensions. 
Namely, we use the fact that the solutions ip £j i(x±) in 3 + 1 dimensions are obtained from 
the solutions ipe 1 \x±) in 2 + 1 dimensions. Thus, spectra in 3 + 1 dimensions can be 
obtained from the results in 2 + 1 dimensions with the substitution M by M, and the 
relation (38). As a consequence, we obtain an additional interpretation of Figs. 1, 2. In 
particular, Fig.l presents energy lowest levels for particles with spin s — 1, and Fig. 2 
presents energy lowest levels for particles with spin s = — 1. 



4 Solenoid regularization 

One can introduce the AB field as a limiting case of a finite radius solenoid field (the 
regularized AB field). In this way, one can fix the extension parameters. First, the manner 
of doing that in the pure AB field was presented by Hagen [9]. Below, we consider the 
problem in the presence of the uniform magnetic field. To this end we have to study 
solutions of the Dirac equation (1) in the combination of the regularized AB field and the 
uniform magnetic field. 

Let the solenoid have a radius R. We assume that inside the solenoid there is an 
axially symmetric magnetic field B m (r) that creates the flux $ = (Zo + n) $0, ^0 — 27r/e. 
Outside the solenoid (r > R) the field B m (r) vanishes. Thus, e B m (r)rdr — l + fi. 
The function B' m (r) is arbitrary but such that integrals in the functions , b(x) in 
(70) are not divergent. We select the potentials of the field B m (r) in the form 



eAf = d (x) eA 2 = -d (x) ^ , (65) 

Rx Rx 



where 

0(x) 



[ f (x) x'dx', f(x) = R 2 eB in (xR), x = r/R. 
Jo 



The potentials of the uniform magnetic field are 

A = 0, A 1 = A (r) A 2 = -A (r) ^ , A (r) = Br 2 /2 . (66) 

r r 

Outside the solenoid the potentials have the form (3). 

Let us analyze solutions of the Dirac equation in the above defined field. To this end 
we have to solve the equation inside and outside the solenoid and continuously join the 
corresponding solutions. The former Dirac spinors we are going to call the inside solutions, 
whereas the latter ones the outside solutions. 

First, we study the problem in 2 + 1 dimensions. We demand the solutions to be square 
integrable and regular at r — > 0. By the same manner as in the Sect. 2, we can find that 
the inside radial spinors (r) (r < R) obey the equation: 

h m ^S, (r) = eiC, (r) , h m = IT™ + a 3 M , 
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where 



n m = _ j d + _ 



R 



x 



l-k- l -{l-^)+'&(x)+i PR x 2 



a\ p R = 1 R 2 /2 . (67) 



We demand the functions tp^ (r) to be square integrable on the interval (0, R). For lu = 
(\e\ = M) the solutions read, 

= (^i, /-io>i, 

= , J-Jo<o, 



x 



ex'*'' exp 



dxx- 1 (tf(x) + £p fi 5 2 ) L 77 = Z - Z - (1 + <r) /2 , (68) 



where c is an arbitrary constant. For cu ^ we present the spinors in the form 



Ci (0 = 



■02 n (?") 



= [a 3 (e - U in ) + M] [d^x)^ + ic-i^L^-i] , 



where c CT are arbitrary constants. The functions (x) satisfy the equation 



1 d d 



x^-- 2 {ri + "d\x)+ ip R x 2 ) 2 + uR 2 -a (f (x) + 2£p R 



x dx dx x 



(69) 



and must be regular at r = in order to satisfy the square integrability condition for 
are interested in the limiting case R — > 0. For our purposes it is enough to 
use the approximation /5j< 1, cui? 2 <C 1. Dropping terms proportional to i? 2 in (67) and 
(69), we find that solutions of Eq. (69) have the form 



cx M e ab{x \ arj > 0, 

cx -H e <yh{x) J* d ~~2\ v \-l e -2ab(x)^ ar]< Q } 



b(x) = [ dxx _1 tf(x) . 
Jo 

The outside solutions (r > R) obey the equation 

J<? (r) = (r) 



(70) 



(71) 



and must be square integrable on the interval (R,oo). Here h is defined by Eqs. (10), 
(11). The general form of the outside solutions reads: 

C? (r) = [a 3 (e - n) + M] (c^r)^ + ic-x^r)^) , 

(r) = (p) , a = l + p-(l + a)/2, 2\ = uft - £(l + p- (1 - a) /2) (72) 



The solutions ip™f (r) and (r) must be joined continuously at r = i?, 

^,° ut (R) = ^ n (R) 
and obey the normalization relation 

N in , N out 1 



(73) 



«i(0) T <,(r)rdr, 1V# 
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(C?(0) f C?(0^r. 



(74) 



One can treat the AB field as a limiting case of a finite radius solenoid field if 

\imN™=0. (75) 

The joining condition (73) one can realize using the following conditions for the functions 

4>™ a (r) and (r) at r = R, 

4> (R - e) = (f> (R + e) , -^0 (i? - e) = ^0 (i? + e) . (76) 
It is convenient to use in (72) the representation (98) for tpx,a (p)- Then, the functions 

^out ^ read; 

1 — oi \ -\~ a 

<Ka ( r ) = a <x J n CT ,m CT (p) + Mm CT ,n CT (p) , Tl a = A — , m CT = A — , 

*r ■ h u ■ k Vra + ^rq + m^) 

a a = K smn ff ff, = — K smm CT 7r, i\ = — . (77) 

sin (n a — m a ) n 

where n a , m a are real numbers. 

Using (76) one can find the coefficients a CT , b a : for the case Z — Z < 0, 

ax = P^-^V, = p£ + ^ 1)/2 ^, (78) 
a^^-^^cL^.^p^-i, (79) 

whereas for the case I — l > 0, 

0l = ^"-^Vai, 6x = pg f '*- 1)/2+ V^, (so) 

a. 1 = p^ + ^a!_ 11 b^ 1 =p^ 2 c'b'_ 1 , (81) 

where the no n- vanishing coefficients a, 6, a', 6' are not depending on p R , the coefficients 
c, c' are normalizing factors which depend on p R . 

Calculating the normalization factors one obtains at R — > 0, for Z — l < 0, 

ai = const ^ 0, &i = 0, a_i = = 0, / > 1 , 
ai = 0, bi — const ^ 0, a_i = 6_i = 0, / < , 

whereas for / — l > 0, 

ai = bi — 0, a_i = const ^ 0, = 0, I > , 
a i = &i = 0, a_i = 0, &_i = const ^ 0, / < — 1 . 

For / = the value of the coefficients is defined by sgn$. One can verify that the condition 
(75) is satisfied. 

Thus, one obtains that for any sign of B the solutions are expressed via Laguerre 
polynomials (23). Particularly, for I = we find that the solutions ip"^ (r) coincide with 
either tp^ (r) or ip 1 ^ (r) accordingly to sgn (<&), 

./.out M _ J (r) , sgn ($) = +1 
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In Sect. 3 we have found the relation between the extension parameter values and 
solution types in the critical subspace / = (63), (64). Now we are in position to refine 
this relation. Namely, if one introduces the AB field as a field of the finite radius solenoid 
for a zero-radius limit, then the extension parameter 6 is fixed to be 6 = — sgn (<3>) ir/2. 
Besides, this way of the AB field introduction explicitly implies no additional interaction 
in the solenoid core. 

To solve the problem in 3 + 1 dimensions we use the results in 2 + 1 dimensions presented 
above. In the limit R — > solutions in the critical subspace have the form 



out 



N 



i+(> + s m) /m] 9o 



1 + (p 3 + sM) /M g (<p)^ 



'.out 



(83) 



where the functions go{(p), VC? ( r ) are defined in (9) and (82), (23), respectively, 
specify the values of the extension parameters in 3 + 1 dimensions as follows, 



We 



e+i = 0-i = -sgn ($) vr/2 . 



(84) 



5 Summary 

We have studied in detail solutions of the Dirac equation in the magnetic-solenoid field 
in 2 + 1 and 3 + 1 dimensions. In the general case, solutions in 2 + 1 and 3 + 1 dimensions 
are not related in a simple manner. However, it has been demonstrated that solutions in 
3 + 1 dimensions with special spin quantum numbers can be constructed directly on the 
base of solutions in 2 + 1 dimensions. To this end, one has to choose the z-component of 
the polarization pseudovector S 3 as the spin operator in 3 + 1 dimensions. This is a new 
result not only for the magnetic-solenoid field background, but for the pure AB field as 
well. The choice S* 3 as the spin operator was convenient from different points of view. For 
example, solutions with arbitrary momentum p 3 are eigenvectors of the operator S 3 . This 
allows us to separate explicitly spin and coordinate variables in 3 + 1 dimensions. Thus, 
in 3 + 1 dimensions one has to study self- adjoint extensions of the radial Hamiltonian 
only. Moreover, boundary conditions in such a representation do not violate translation 
invariance along the natural direction which is the magnetic-solenoid field direction. The 
self-adjoint extensions of the Dirac Hamiltonian in the magnetic-solenoid field have been 
constructed using von Neumann's theory of deficiency indices. A one-parameter family 
of allowed boundary conditions in 2 + 1 dimensions and a two-parameter family in 3 + 1 
dimensions have been constructed. By that the complete orthonormal sets of solutions 
have been found. The energy spectra dependent on the extension parameter have been 
defined for the different self-adjoint extensions. Besides, for the first time solutions of the 
Dirac equation in the regularized magnetic-solenoid field have been described in detail. 
We considered an arbitrary magnetic field distribution inside a finite-radius solenoid. It 
was shown that the extension parameters G = — sgn($)7r/2 in 2 + 1 dimensions and 
+ i = 0_i = — sgn($)7r/2 in 3 + 1 dimensions correspond to the limiting case R — > of 
the regularized magnetic-solenoid field. 
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A Appendix 

l.The Laguerre function I n , m (x) is defined by the relation 

InM = \/ w ( !^1rn P i~ X/2 \ " ( "" m)/2 ^-^ n-m + l;x). (85) 
y r (1 + m) r (l + n — m) 

Here <3> (a, b; x) is the confluent hypergeometric function in a standard definition (see [31], 
9.210). Let m be a non-negative integer number; then the Laguerre function is related to 
Laguerre polynomials L^(x) ([31], 8.970, 8.972.1) by the equation 



Wv»(z) = y r(m ^ Q + 1) ^ / V/ 2 L^(x), (86) 
1 d m 

L°(x) = —e x x- a —e~ x x m+a . (87) 

mK 1 m\ dx m v 1 

Using well-known properties of the confluent hypergeometric function ( [31], 9.212; 9.213; 
9.216), one can easily get the following relations for the Laguerre functions 

2y/x(n + l)I n+ljTn (x) = (n-m + x)I n , m (x) - 2xl' nm (x) , (88) 

2y/x(m+ l)I njm+1 (x) = (n-m- x)I n , m (x) + 2xl' nm (x) , (89) 

2y/xnl n - 1 , m {x) = (n-m + x)I n , m (x) + 2xl' nm (x) , (90) 

2^/xmI n ^i(x) = (n-m- x)I njTn (x) - 2xl' nm (x). (91) 

Using properties of the confluent hypergeometric function, one can get a representation 



,r(l+n) exp(a;/2) n-m , , _ , . . 

I n ,m(x) = t r i , \ rn ^ 1 ' s x— $ 1 +n ,l + n-m;-x), 92 
w r(l + m) r(l + n — m) 

and a relation ([31], 9.214) 

In,rn( x ) = (~ 1 T~ mi m,n (x) , n - m integer . (93) 
The functions / a+mjm (:r) obey the ortho normality relation 

/•oo 

/ Ia+n,n (x^) Ia+rn;rn (x) dx &m,n i (^4) 
./O 
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which follows from the corresponding properties of the Laguerre polynomials 
7.414.3). The set of the Laguerre functions 

I a +m, m (x), m = 0, 1, 2... , a > -1 

is complete in the space of square integrable functions on the half-line (x > 0), 

oo 

^] Ia+m,m( x )Ia+m,m(y) = $ ( x ~~ U) ■ 

2. The function if)\ >a (x) is even with respect to index a, 

^A,a (x) = I^X-a (x) . 

It can be expressed via the confluent hypergeometric functions 

T (—a) x% ( 1 + a 



ip\, a (x) = e 2 

, r(q)x-f 
+ r(^-A)"V 2 

or, using (85), via the Laguerre functions 



r(^-A) 



A, 1 + a; x 



<]> | 1 _ Q - A, 1 - a; x 



ip\,a (x) = — ^ + n ) r ^ + m ) (sinn7r/ ni „ t (x) - smmirl m:n (x)) , 
sin (n — m)n 

1 — at 1 + a 
a = n — m, 2A = 1 + n + m, n = \ — , m = A — . 

There are the following relations of the functions ip\ >a (x) , 



[31], 



(95) 



(96) 



(97) 



(98) 



/ / \ /— , / 1 + a — 2A , . . 

l/j\,a{x) = M-i^-l( I ) + 2 ^-l,aW, 

^A.afa) = V^A-i.a+1 0*0 + ^ ^ A-l.a (x) , 

2^ A , a (x) = (2A - 1 - x) *Px, a (x) + l - (2A - 1 - a) (2A - 1 + a) ^ A _ 1)Q (x) , 

2 ^A,a (x) = (a-x) iJ X ,a (x) + (2A - 1 - a) V^A-I.a+l (&) 
= (x - 2A - 1) ^A,a - 2^ A +l,a ■ 

As a consequence of these properties we get 

2A- 1 + a 



A a ^ X ,a (x) = ^ ' "A -j.a-l (x) , A+^A-i.a-l (&) = (&) , 

x + a ^ ^j— d ^ + x + a — 1 d 
2y / X cix' a 2y / x dx 



Using well-known asymptotics of the Whittaker function ([31], 9.227), we have 

(x) ~ X A_ 5 e -f , X ^ OO; ^A,a (x) - 



(i^-A) 



(99) 



(100) 



r ( I Cej) \a\ , 

a;"-, a + 0, x ~ 0. (101) 
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The function i/j\ tOC (x) is correctly defined and infinitely differentiable for < x < oo and 
for any complex A, a. In this respect one can mention that the Laguerre function are not 
defined for negative integer n, m. In particular cases, when one of the numbers n, m is 
non-negative and integer, the function ip\ t0l (x) coincides (up to a constant factor) with 
the Laguerre function. 

According to (101), the functions ip\ t a (x) are square integrable on the interval < x < 
oo whenever \a\ < 1. It is not true for \a\ > 1. The corresponding integrals at a ^ can 
be calculated as following ([31], 7.611), 



oo 

/ 



oo 

/ 



1p\,a (x) lj>x, a (x) dx 



1 - a - 2A' 



\lp X ,a (x) | 2 dx = 



IT 



(A' — A) sinew 



1 + a - 2A 



l + a-2X' 



\a\ < 1, 



1 - a - 2A N 



|a| < 1, 



- -i 

(102) 
(103) 



Here tp(x) is the logarithmic derivative of the T-function ( [31], 8.360). In the general 
case, the functions ip\ a (x) and ipy a (x), A' ^ A, are not orthogonal, as it follows from 
(102). 



B Appendix 

Below we present modifications of some above formulas for the case B < 0. 

1. The spectrum of u corresponding to the functions m ,« i(T (r) is 

f2 7 (m-I + l-^I-(l + a)/2<0 

\ 2 7 (m + (l-<7)/2), Z-(l + <r)/2>0 ' 1 j 

and the spectrum of c<j corresponding to the functions </>^ CT (r) is 

y 27m, (7 — 1 

These expressions are modifications of Eqs. (18), (19) for the case B < 0. 

2. Consider the spinors (r) satisfying (20). In the case ui = they read 

«o^(r) = ( ^' l(r) ),/>!; ttf(r) = ( ^' l(r) ) , Z = 0. (106) 
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In the case to ^ they are 

= ( ^,,r) ) ■ ' < "I , - = " ( + 1 " /<) , 

»-u, ± M = ( ),l>i,« = 2 T (ro + i), 



r 



<±M = ( ^^V) ) ' ' = • 10 = 27 (m + 1 " ' (107) 

These expressions are modifications of Eqs. (21), (22) for the case B < 0. 

3. In the case u = only positive energy solutions (particles) of Eq. (10) are possible. 
They coincide with the corresponding spinors u up to a normalization constant: 

+ ?Mr) = Nu 0il (r), I > 1; +^(r) = Nu^r), I = 0. (108) 

Thus, particles have the rest energy level, and the antiparticle states spectrum begins 
from = -a/M 2 + 27. 

4. Relations for the irregular spinors u Ui(T (r), similar to ones (27), for the case B < 
have the form 

nu w i(r) = i^/2r/u w ,i(r), Uu W)1 {r) = -i-^=« w i(r) . (109) 
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